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SECTION  I 
INTRODUCTION 


The  purpose  of  the  research  presented  in  this  document  is  to  refine  our  ability  to  ! 

discriminate  between  nuclear  explosions  and  earthquakes  based  on  the  depth  of  an  event  I 

determined  from  local  and  regional  phases.  If  the  depths  of  seismic  events  can  be  deter¬ 
mined  accurately,  many  events  can  be  eliminated  as  potential  nuclear  explosions  based  ! 

on  the  practical  limit  of  nuclear  explosion  burial  depth.  This  discriminant  has  been  used 
for  many  years  at  teleseismic  distances.  The  higher  frequency  data  obtained  at  local  and 
regional  distances  improves  the  depth  resolution  of  the  method  and  thus  the  minimum  j 

depth  at  which  the  discriminant  may  be  applied.  In  addition,  smaller  events  not 
recorded  at  teleseismic  distances  may  be  seen  at  regional  and/or  local  distances. 

Our  approach  to  this  problem  has  been  to  improve  phase  identification  and  picking 
abilities  of  depth  phases  pP  and  sP  through  filtering  and  three-component  polarization 
methods.  These  methods  are  applied  to  data  from  a  set  of  earthquakes  located  in  the  I 

northeastern  United  States  and  adjacent  Canada.  The  area  of  study  was  chosen  for 
several  reasons: 

(1)  the  large  majority  of  events  in  this  area  occur  at  depths  less  than  15  km  so  that  a 
data  set  of  shallow  events  is  readily  available, 

(2)  the  geology  and  tectonic  setting  of  the  area  is  well  known  and  similar  to  that  found 
in  the  Soviet  Union, 

(3)  the  depths  of  the  events  in  the  area  of  study  are  fairly  well  known  from  network 
studies  and  aftershock  surveys,  and 

(4)  there  are  broad-band  three-component  digital  seismic  stations  in  the  immediate  I 

vicinity  (RSNY  and  SRNY). 

The  most  difficult  problem  associated  with  depth  determination  at  local  and 
regional  distances  using  depth  phases  is  identification.  As  the  depth  phases  are  usually 
buried  in  the  coda,  data  processing  plays  a  key  role  in  identifying  and  picking  them. 

Band-pass  filtering  and  three-component  polarization  state  filtering  methods  are  used  on  ! 

the  data  in  an  attempt  to  isolate  the  depth  phases  and  adaptive  polarization  analysis  is 
used  to  help  identify  them. 

Once  the  depth  phases  have  been  correctly  identified  and  picked,  the  depth  can  be 
determined  accurately  if  the  source  region  velocity  structure  is  known.  However, 
misidentification  of  phases  and  inaccurate  time  picks  can  lead  to  severe  errors  in  the 
depth  estimates.  To  assess  the  potential  depth  errors  resulting  from  these  problems,  syn¬ 
thetic  seismogram  depth  sections  have  been  generated  for  each  of  the  events  to  compare 
to  the  identification  and  time  picks  made  from  the  data.  Of  course,  the  depth  of  the 
event  is  not  the  only  factor  contributing  to  the  character  of  the  depth  section;  the 
velocity-depth  structure  used  and  the  focal  mechanism  are  also  very  important  to  the 
overall  character  of  the  section. 
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This  report  covers  the  first  year’s  effort  of  a  two  year  research  project.  In  the  sections 
that  follow,  we  briefly  describe  the  analysis  methods  used  on  the  data  to  isolate  depth 
phases,  review  the  set  of  earthquakes  that  we  have  selected  for  study,  present  the  syn¬ 
thetic  seismograms  for  these  events,  and  show  a  detailed  analysis  of  the  January  19,  New 
Hampshire  earthquake. 
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SECTION  II 

POLARIZATION  ANALYSIS  METHODS 


Polarization  analysis  of  a  three-component  wavefield  is  used  to  determine  the  vector 
particle  motion  characteristics  such  as  the  degree  of  polarization,  direction  of  propaga¬ 
tion  and  ellipticity.  This  information  can  then  be  used  to  design  filters  that  emphasize 
or  eliminate  waves  having  specific  polarization  properties.  In  seismology,  this  method 
has  been  applied  to  analyze  and  filter  waves  coming  from  certain  directions,  or  having 
rectilinear  or  elliptical  particle  motion,  and  to  discriminate  against  isotropic  noise  (Sut¬ 
ton  and  Pomeroy,  1963;  White,  1964;  Archambeau  et  al.,  1965)  . 

Polarized  arrivals  can  be  enhanced,  their  wave  type  identified,  and  their  direction  of 
propagation,  as  well  as  their  phase  velocity  for  body  waves  and  ellipticity  for  Rayleigh 
waves,  determined.  This  information  can  then  be  used  for  phase  identification,  epicen- 
tral  location,  source  depth  determination,  earth  structure  studies,  ...etc.  Data  filtered  for 
polarized  signal  are  more  amenable  to  comparison  with  synthetic  seismograms,  providing 
a  better  determination  of  relevant  parameters. 

The  polarization  state  method  was  proposed  by  Samson  (1977)  for  the  analysis  and 
filtering  of  seismograms.  It  makes  use  of  the  spectral  matrix,  which  is  the  frequency 
domain  equivalent  of  the  cross-correlation  matrix  of  the  data,  smoothed  over  a  frequency 
window.  Various  spectral  estimators  can  be  obtained  from  the  spectral  matrix,  such  as 
the  degree  of  polarization,  the  degree  of  rectilinear  polarization,  a  detector  for  signal  in  a 
specific  "pure  state"  (i.e.  completely  polarized)  or  with  specific  polarization  characteris¬ 
tics.  Such  estimators  can  then  be  used  to  filter  the  data  in  order  to  discriminate  against 
unpolarized  noise  or  to  enhance  signal  in  some  specific  polarization  state  (Samson  and 
Olson,  1981).  Since  the  polarization  characteristics  of  seismograms  continuously  change 
with  time,  such  filters  have  to  be  data-adaptive:  a  sliding  time  window  of  data  is 
analyzed  and  a  filter  is  specifically  designed  for  that  window.  Polarization  properties  are 
determined  and  used  not  only  as  a  function  of  time,  but  also  as  a  function  of  frequency, 
since  the  processing  is  done  in  the  frequency  domain.  This  is  an  advantage  that  methods 
which  operate  purely  in  the  time  domain  do  not  have. 

To  enhance  depth  phases  pP  and  sP  we  take  advantage  of  the  three-component 
polarization  state  filter’s  ability  to  pass  frequencies  that  contain  rectilinearily  polarized 
motion  and  reject  frequencies  that  do  not. 

After  the  data  have  been  state  filtered  for  rectilinear  motion,  they  are  processed 
with  an  adaptive  polarization  method.  The  horizontal  components  are  adaptively 
rotated  to  the  time  varying  radial  and  transverse  directions  and  the  time  variable 
azimuth  and  apparent  angle  of  incidence  are  given.  This  process  helps  in  the 
identification  of  phases.  Also,  the  synthetic  seismogram  components  are  the  radial  and 
vertical  so  rotating  the  horizontal  data  to  the  adaptive  radial  and  transverse  components 
facilitates  comparison. 
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The  polarization  state  filter  and  the  adaptive  polarization  process  form  the 
the  data  processing  for  this  project. 
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SECTION  III 
NORTHEAST  EVENTS 


Eight  local  and  regional  earthquakes  and  one  quarry  blast  were  selected  for  com¬ 
parison  with  synthetic  seismograms  discussed  below.  The  events  were  chosen  for  their 
size,  location,  and  availability  of  data.  A  list  of  these  events  is  given  in  Table  1  along 
with  location,  depth,  magnitude,  strike,  dip,  and  rake  of  the  focal  mechanism,  the  region 
(either  New  England  or  Grenville),  and  references  for  the  focal  parameters.  The  event 
and  station  locations  are  plotted  in  Figure  1,  showing  their  spatial  distribution.  Many  of 
the  earthquakes  in  Table  1  have  been  studied  by  others  and,  for  some,  the  focal  depths 
and  focal  mechanisms  are  well  determined  and  will  enable  us  to  test  the  methods  we  are 
developing. 

In  addition  to  reviewing  the  literature,  we  examined  phase  data  from  local  seismic 
networks  originally  used  to  invert  for  hypocentral  location.  In  three  cases  the  estimated 
error  on  focal  depth  was  so  large  that  to  use  the  reported  depth  would  be  misleading. 
For  these  we  indicate  (Table  l)  an  unknown  depth. 

Figures  2A  through  21,  all  in  essentially  the  same  format,  show  the  three-component 
seismograms  for  events  in  Figure  1.  Most  were  recorded  by  the  Regional  Seismic  Test 
Network  (RSTN)  station,  RSNY,  in  northern  New  York  and  several  of  the  recent  events 
have  been  recorded  by  SRNY,  a  very-broad-band  seismic  station  operated  by  Rondout 
Associates  in  Stone  Ridge,  New  York.  Just  the  P- wave  portion  of  the  seismograms  is 
shown  in  these  figures  since  depth  phases  are  our  primary  interest.  As  mentioned  in  the 
introduction,  one  of  the  keys  to  understanding  the  data  and  being  able  to  compare  them 
to  synthetic  seismograms  will  be  the  extraction  of  discrete  phase  arrivals  from  the  P- 
coda  waves.  Along  with  the  unfiltered  or  high-pass  filtered  data  (the  top  three  traces  in 
figures  2A  through  21),  we  present  three-component  state  filtering  of  the  data  (the  bot¬ 
tom  three  traces).  These  are  filtered  to  pass  rectilinear  particle  motion  as  a  function  of 
both  frequency  and  time  (as  described  in  analysis  methods  above).  In  all  cases, 
polarization-state  filtering  serves  to  enhance  some  portions  of  the  seismograms  relative  to 
other  portions.  We  are  currently  investigating  ways  of  identifying  the  phases  in  these 
"cleaned  up"  seismograms.  An  example,  the  January  19,  1982  New  Hampshire  earth¬ 
quake,  is  presented  following  the  discussion  of  synthetic  seismograms. 
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New  Horpshire 
19  January,  1982  RSNY 

a=  267  Wn  « 

Figure  2A.  Seismograms  recorded  at  RSNY  of  the  January  19,  1982,  New  Hampshire 
earthquake.  The  top  three  traces  are  the  vertical,  north-south, 
and  east-west  components,  unfiltered.  The  bottom  three  traces 
are  the  three  components  after  state-filtering  for  rectilinear 
particle  motion.  The  rest  of  the  events  shown  in  Figure  1  are  4 
similarly  depicted  in  Figures  2B  through  21. 
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Ontario 

11  October,  1985  RSNY 
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Figure  2E.  11  October  1983  Ontario  earthquake. 


Adirondack 

23  October,  1984  RSNY 
a=  115  m 


Figure  2F.  23  October  1984  Adirondack  earthquake. 
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Figure  2G.  19  October  1985  Ardsley  earthquake. 
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SECTION  IV 

SYNTHETIC  SEISMOGRAMS 


The  synthetic  seismograms  were  computed  using  the  locked  mode  method  (Harvey, 
1981).  Frequencies  from  zero  to  5  Hz  are  computed.  Regional  geology  in  the  Northeast 
required  that  two  different  velocity/ attenuation  models  be  used  for  the  computation;  a 
New  England  model  for  the  events  in  New  England  and  southern  New  York,  and  a  Gren¬ 
ville  model  for  events  within  the  Grenville  province  (Figure  1,  Table  2).  The  New  Eng¬ 
land  velocity  model  was  derived  by  Taylor  et  al.,  (1980)  using  regional  travel  times  of  P 
and  S  waves  recorded  across  the  Northeastern  United  States  Seismic  Network  operated 
by  several  academic,  governmental  and  private  institutions  in  the  northeast.  The  Gren¬ 
ville  model  was  from  unpublished  refraction  results.  The  frequency  and  depth  dependent 
Q  model  for  the  New  England  model  was  adapted  from  work  by  Mitchell  (1981)  and  a 
simple  frequency  dependent  Q  was  used  for  the  Grenville  model.  Although  there  are 
significant  differences  in  the  Q  models,  there  is  little  effect  on  the  synthetics  because  the 
frequencies  and  ranges  at  which  they  are  being  computed  are  relatively  small. 

Synthetics  were  computed  for  each  event  at  several  depths  with  the  reported  focal 
mechanism  listed  in  Table  1.  The  instrument  response  was  convolved  with  the  synthet¬ 
ics  for  comparison  to  the  data. 

Station  RSNY  lies  within  the  Grenville  structure  close  to  the  boundary  with  the 
New  England  Structure.  While  the  travel  paths  to  RSNY  from  New  England  events 
must  pass  through  both  structures,  the  majority  of  the  path  is  within  the  New  England 
structure  and  so  the  New  England  structure  is  a  good  approximation  to  the  true  path. 

In  Figures  3 A  through  3H  we  show  the  synthetic  depth  sections  for  the  events  in 
Table  1.  Only  the  first  20  or  so  seconds  are  shown  as  we  wish  to  emphasize  the  charac¬ 
teristics  of  the  depth  phases  with  focal  mechanism,  depth,  range,  and  velocity  model. 
Some  of  the  phases  are  identified  in  the  synthetic  sections  with  labeled  lines.  Reflections 
off  of  the  13  km  discontinuity  in  the  New  England  model  and  the  4  km  discontinuity  in 
the  Grenville  model  are  indicated  with  an  "i”  (e.g.  PiP ).  Mantle  reflections  are  indi¬ 
cated  by  an  "m".  The  depth  phases  are  characterized  by  progressively  later  arrivals  with 
increased  depth  and  the  primary  phases  by  either  very  little  change  in  arrival  time  with 
depth  for  the  direct  arrivals  or  progressively  earlier  arrivals  with  increased  depth.  Notice 
that  in  some  of  the  sections,  the  depth  phases  dominate,  while  in  others,  the  primary 
phases  dominate.  This  is  a  function  of  the  theoretical  radiation  pattern  predicted  by  the 
input  focal  mechanism.  Because  the  input  velocity  structure  will  influence  the  take-off 
angle  of  each  phase  and  hence  its  position  on  the  focal  sphere,  both  focal  mechanism  and 
velocity  model  are  contributing  to  the  relative  amplitudes  of  primary  and  depth  phases. 
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Figure  3B.  Synthetic  seismogram  depth  section.  Source®  31  August 
1982  Earthquake.  Strike®  173,  Dip®  60,  Rake®  110. 
Receiver  response®  RSNY  (sp) .  Azimuth®  330. 
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Maine  Synthetics 
29  May,  1983 
New  England  Velocity  Model 
a=  3 22  KJn 

Figure  3C .  Synthetic  seismogram  depth  section.  Source=  29  May,  1983 
Earthquake.  Strike*  355,  Dip  =  52,  Rake=  78.  Receiver 
response*  RS'JY  (sp)  .  Azimuth*  269. 
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Grenville  Velocity  Model 
a=  69  Wn 


Figure  3D.  Synthetic  Seismogram  depth  section.  Source=  7  October 
1983  Earthquake.  Strike=  173,  Dip=  60,  Rake=  110. 
Receiver  response=  RSNY  (sp) .  Azimuth=  342. 
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Adirondack  Synthetics 
23  October,  1984 
Grenville  Velxity  Model 
a=  113  Wn 


Synthetic  Seismogram  depth  section.  Source”  23  October 
1984  Earthquake.  Strike*  173,  Dip*  60,  Rake*  110. 
Receiver  response*  RSNY  (sp) .  Azimuth*  335. 
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Ards ley  Synthetics 
19  October  1985 
New  England  Velocity  Model 
a=  100  Wn 

Figure  3G.  Synthetic  Seismogram  depth  section.  Source®  19  October 
1985  Earthquake.  Strike*  022,  Dip*  90,  Rake*  180. 

Receiver  response-  SRNY .  Azimuth*  344.  Radial  components 
shown . 
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New  England  Velocity  Model 
a=  120  Km 


Figure  3H.  Synthetic  Seismogram  depth  section.  Source  at  zero  depth= 

2  October,  1985  Explosion.  Source  at  other  depths®  30  October 
1985  Earthquake.  Strike®  18,  Dip=  60,  Rake®  156.  Receiver 
response®  SRNY.  Azimuth®  180  from  explosion  and  181  from 
Earthquake.  The  blast  and  earthquake  epicenters  are  only 

3  km  apart. 
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TABLE  2.  VELOCITY/Q  MODELS 
New  England  Model  (Taylor  et  al. ,  1980) 
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SECTION  V 

ANALYSIS  OF  JANUARY  19,  1982  NEW  HAMPSHIRE  EARTHQUAKE 

The  success  of  determining  focal  depth  using  depth  phases  and  comparing  them  to 
synthetics  depends  largely  on: 

(1)  good  approximations  of  crustal  velocity  structure 

(2)  recovery  and  correct  identification  of  the  phases  in  the  data. 

In  this  section  we  use  one  earthquake  to  illustrate  how  the  synthetic  seismograms  are 
affected  by  2  different  velocity  models  and  focal  mechanisms  as  well  as  to  illustrate  our 
approach  to  identifying  depth  phases. 

Influence  of  some  Input  Parameters  on  Synthetics 

The  earthquake  chosen  is  the  Gaza,  N.H.,  magnitude  4.5  event  of  January  19,  1982. 
Since  the  source  is  in  the  Appalachian  Province  and  the  receiver  (RSNY)  is  in  the  Gren¬ 
ville  Province  (see  Figure  1)  we  calculated  separate  synthetics  using  the  two  velocity 
models  (Table  2)  to  test  which  is  more  appropriate. 

Figure  4  shows  radial  and  vertical  components  for  the  initial  15  seconds  of  the  syn¬ 
thetic  seismograms  calculated  for  two  velocity  models.  RSNY  is  the  theoretical  receiver, 
a  distance  of  267  km  and  an  azimuth  of  295°  from  the  New  Hampshire  source.  The  velo¬ 
city  models  and  Moho  reflections  are  shown  schematically  for  comparison.  The  focal 
mechanism  is  given  by:  strike  =  280°,  dip  =  75°  and  rake  =  -11°.  The  input  focal 
depths  are  comparable. 

As  expected,  the  prominent  phases  and  their  arrival  times  are  different  for  the  two 
different  velocity  models.  The  Pn  phase  is  nodal  for  the.  radiation  pattern,  despite  small 
differences  in  take-off  angles  for  the  two  models.  It  is  simply  indicated  at  the  appropri¬ 
ate  arrival  time.  For  the  New  England  model,  PmP  is  small  relative  to  sPmP.  They 
both  arrive  before  the  intra-crustal  reflections  PiP  and  sPiP .  For  depths  greater  than  10 
km,  however,  PiP  arrives  before  sPmP  (see  figure  3A).  PiP  and  sPiP  are  sharp  arrivals 
because  of  the  large  velocity  contrast  at  15  km  in  the  crustal  model.  Vp  increases  from 
6.1  km/sec  to  7.0  km/sec  at  this  boundary. 

For  the  Grenville  model,  Pg  precedes  the  Moho  reflections  (for  a  source  depth  of  7 
km).  The  amplitude  of  sPmP  relative  to  PmP  is  smaller  for  the  Grenville  model  than 
for  the  New  England  model  because  the  upgoing  sPmP  is  crossing  a  fairly  sharp  internal 
discontinuity  near  the  source  while  the  downgoing  PmP  only  crosses  this  boundary  once, 
near  the  receiver.  In  the  case  of  New  England  structure,  the  sharp  discontinuity  is  below 
the  source,  instead  of  above  it.  In  addition,  PmP  is  close  to  a  nodal  plane  and  small 
differences  in  take-off  angle,  resulting  from  differences  in  the  velocity  structures,  may  be 
contributing  to  the  difference  in  the  amplitudes  of  PmP. 
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January  19 ,  1982,  a=  267  Km,  Azimuth=  295. 

Figure  4.  Two  velocity  models  and  resulting  synthetic  seismograms. 


A*  New  England  Model,  source  depth*  6  Km.  Source-to-receiver  Moho,  reflections 
shown  schematically.  15  seconds  of  radial  and  vertical  component  synthetic^ 
with  phases  identified.  PiP  and  sPiP  are  internal  reflections  at  the 
15  Km  boundary. 

B*  Grenville  model,  source  depth*  7  Km.  Source-to-receiver  Moho  reflections 
shown  schematically.  15  seconds  of  radial  and  vertical  component  synthetics 
with  phases  identified. 
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Even  though  the  two  velocity  models  produce  strikingly  different  synthetic  seismo¬ 
grams,  it  is  important  to  note  that  for  the  depths  and  models  shown,  the  time  between 
the  arrivals  of  the  downgoing  and  upgoing  Moho  reflections  is  the  same.  Thus,  if  one 
can  independently  identify  these  two  phases  in  the  actual  data,  a  match  (ignoring  abso¬ 
lute  travel  times)  to  the  Grenville  synthetic  would  suggest  a  depth  of  7  km  and  a  match 
to  the  New  England  synthetic  would  suggest  a  depth  of  6  km.  Given  this  time 
difference  between  the  two  arrivals,  the  choice  of  velocity  model  has  had  only  a  small 
effect  on  depth  determination.  If,  however,  one  uses  the  computed  phase  travel  times  to 
help  identify  depth  phases  in  the  real  data,  the  variation  in  travel  time  with  velocity 
model  could  lead  to  misinterpretations  of  the  real  data. 

Next,  in  Figure  5  we  show  15  seconds  of  synthetics  calculated  for  2  different  focal 
mechanisms.  The  mechanism  shown  at  the  bottom  was  determined  from  P- wave  first 
motion  data,  mainly  from  the  Northeast  United  States  Seismic  Network  (Pulli  et  al., 
1983).  The  mechanism  at  the  top  of  the  figure  differs  mainly  in  the  quadrants  of  dilata¬ 
tion  and  compression.  The  azimuth  to  station  RSNY  is  very  close  to  a  nodal  plane  in 
both  instances,  so  there  is  not  a  great  difference  in  radiation  pattern.  Thus,  the  same 
initial  phases  appear  on  both  sets  of  seismograms,  but  the  amplitudes  of  these  phases  are 
different.  If  theoretical  radiation  patterns  can  be  used  to  predict  actual  phases  in  the 
data,  then  it  will  be  helpful  to  have  independently  determined  fault  plane  solutions. 


Comparison  of  Data  to  Synthetics 

It  would  be  very  difficult  to  compare  the  raw  data  recorded  at  RSNY  to  the  syn¬ 
thetic  seismograms  (compare  the  top  of  figure  2A  with  the  synthetics  in  figure  3A).  With 
the  polarization  filtering,  however,  we  have  been  able  to  recover  polarized  phases  (bot¬ 
tom  of  figure  2A)  that  were  all  but  lost  in  the  P-coda  of  the  unfiltered  seismograms. 
After  polarization  filtering,  we  low-pass  filter  the  data  to  remove  frequencies  greater  than 
5  Hz,  which  are  not  modeled  by  the  synthetics. 

The  following  figures  show  comparisons  of  the  vertical-component  RSNY  processed 
data  to  the  vertical-component  synthetic  depth  section  calculated  for  the  New  England 
model.  Twenty  seconds  are  shown.  We  began  by  inserting  the  data  at  a  depth  of  about 
9  km  (Figure  6).  There  are  several  arrivals  in  the  data  that  match  predicted  arrivals. 
These  are,  in  order,  PmP,  sPmP,  the  dual  arrivals  PiP  and  Pg,  and  finally  an  upgoing 
phase  arriving  toward  the  end  of  this  time  window.  Clearly,  the  match  is  not  perfect. 
Signal  arrives  in  the  data  later  than  the  appropriate  time  for  sPiP,  and,  subsequently 
until  the  late  phase  (in  the  data,  about  14  seconds  after  the  initial  P),  there  is  not  a 
good  correlation  between  data  and  synthetics.  A  comparison  of  travel  times  for  the  first 
arriving  phase  in  the  data  and  synthetics  suggests  that  they  may  not  be  the  same  phase. 
The  travel  time  of  synthetic  PmP  at  a  source  depth  of  9  km  (with  the  given  input  velo¬ 
city  model)  is  40.40  seconds,  whereas  the  travel  time  of  the  first  phase  in  the  data  is 
38.59  seconds.  Allowing  for  uncertainties  in  the  actual  origin  time  of  the  earthquake  and 
for  uncertainties  in  crustal  structure,  they  could  be  the  same  phase.  If,  however,  we 
assume  that  the  origin  time  and  velocity  model  are  appropriate,  we  can  match  several 
early  P-phases  by  aligning  the  data  at  a  depth  of  approximately  4  1/2  km  as  illustrated 
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New  Hampshire  earthquake,  synthetic  with  data.  Same  as  Figure  )A 
except  the  real  Data  are  included  in  the  depth  section  at  a  depth 
of  about  9  Km.  See  text  for  discussion. 
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Figure  7, 


New  Hampshire  earthquake,  synthetics  with  data.  Like  Figure  6  except 
notice  that  the  real  data  are  plotted  at  a  depth  of  about  4.5  Km. 

See  text  for  discussion. 
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in  Figure  7.  Since  the  real-data  travel  time  is  more  appropriate  for  Pn,  the  first  arrival 
is  now  aligned  with  the  predicted  Pn  arrival  time  (compare  Figures  6  and  7).  For  this 
shallower  depth,  the  data  match  predicted  arrival  times  for  Pn,PmP  ,sPmP  and  for  the 
late  up-going  phase.  Again,  the  correlation  is  not  perfect. 

At  this  point,  it  is  clear  that  a  great  deal  depends  on  the  correct  identification  of  the 
arrivals  in  the  data.  In  addition  to  travel  times,  the  azimuth  and  apparent  angles  of 
incidence  can  be  used  to  help  identify  phases.  In  Figure  8  we  show  the  results  of  an 
adaptive  polarization  analysis  of  the  first  20  seconds  of  the  same  state-filtered  seismo¬ 
grams  that  we  have  been  comparing  to  the  synthetics.  The  first  three  traces  in  the  figure 
are  the  vertical,  and  the  adaptive  radial,  and  transverse  horizontal  components.  The 
direction  of  maximum  signal  strength  as  a  function  of  time  represented  by  the  azimuth 
and  apparent  angle  of  incidence,  is  indicated  (with  error  estimates)  below  the  RZ  pro¬ 
duct  trace  in  Figure  8.  The  positive  values  for  the  product  of  the  radial  and  vertical 
components  indicate  that  all  the  phases  are  arriving  at  the  receiver  as  P- waves.  An 
adaptive  polarization  analysis  of  the  synthetic  seismograms  indicated  the  same  thing; 
thus,  at  the  very  least,  we  know  that  the  last  conversions  were  to  P-waves  for  both  the 
real  and  synthetic  data. 

Interestingly,  we  noted  that  the  last  arrival,  which  had  matched  a  predicted  phase 
in  the  synthetics  at  both  the  4.5  and  9  km  depth  positions,  is  probably  not  even  modeled 
by  the  synthetics.  Though  highly  polarized  and  hence  a  prominent  phase,  it  is  fourteen 
degrees  off  the  correct  azimuth  from  the  source.  Unlike  the  two  other  out-of-line 
azimuths  (see  Figure  8),  the  error  estimate  does  not  allow  enough  uncertainty  to  put  it 
in  line.  Furthermore,  this  late  phase  is  approaching  the  receiver  at  a  shallower  angle 
than  any  other  arrival  (see  the  apparent  angles  of  incidence),  yet  the  late  phase  in  the 
synthetic  seismograms  was  approaching  the  receiver  at  the  steepest  angle  in  the  same  20 
second  time  window.  A  steep,  late  arrival  would  be  consistent  with  a  multiply-reflected 
phase. 

The  phase  immediately  proceeding  the  last  phase  (by  almost  1  second)  is  more  con¬ 
sistent  with  the  synthetics  because  that  phase  is  approaching  at  a  steeper  angle.  Using 
that  approach,  a  good  match  is  obtained  by  inserting  the  data  at  a  depth  of  6  km  (Fig¬ 
ure  9).  Here,  the  first  arrival  is  a  predicted  PmP  and  there  are  phases  matching  the 
predicted  sPmP  and  the  dual  arrival  of  Pg  and  PiP,  as  well  as  the  late  multiply- 
reflected  phase.  In  addition,  the  relative  amplitudes  match  the  synthetics  fairly  well. 

We  have  illustrated  some  of  the  difficulties  in  identifying  and  correlating  depth 
phases  from  a  single  set  of  3-component  data.  Evidently,  this  particular  earthquake  is 
difficult  to  pin  down  as  the  range  of  estimates  for  its  focal  depth  is  3-11  km,  using  a 
variety  of  approaches  (see  Table  3). 
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Figure  8.  Adaptive  polarization  of  state-filtered  seismogram  recorded  at 

RSNY.  The  values  of  azimuth  and  apparent  angle  of  incidence  are 
given  with  estimated  errors. 


January  19,  1982 
Gaza,  New  Hampshire 
Vertical  Component  Synthetics 
New  England  Velocity  Model 
A=  267  Km 

Figure  9.  New  Hampshire  earthquake,  synthetics  with  data.  Like  Figures  6 
and  7  except  the  real  data  are  plotted  at  a  depth  of  6  Km  in  i.he 
depth  section.  See  text  for  discussion. 
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TABLE  3.  DEPTH  ESTIMATES  FOR  THE 
GAZA,  NEW  HAMPSHIRE,  EARTHQUAKE  JANUARY  19,  1982 


Depth  Estimate  (km) 

Method 

Reference 

3 

local  network 

Pulli  et  al. ,  1983 

3.5 

teleseismic  body 
wave  modeling 

Pulli  et  al. ,  1983 

9 

relative  relocations 
using  aftershocks* 

Brown  and  Ebel,  1985 

4-11 

regional,  surface 
wave  modeling 

Hermann,  pers. 
Communication 

*  Depth  range  of  larger  aftershocks  recorded  by  portable  instruments  is 
2.7-4. 7km. 
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SECTION  VI 
FUTURE  WORK 

The  work  presented  in  this  report  represents  the  first  years’  research  of  a  two  year 
project.  The  major  portion  of  the  seismogram  synthesis  has  been  completed  and  in  the 
months  ahead  only  minor  changes  to  the  focal  mechanisms  may  be  necessary  for  the 
depth  section. 

We  will  continue  the  detailed  comparisons  of  state-filtered  and  adaptively  polarized 
data  to  the  synthetic  depth  sections  presented  here.  In  addition,  preliminary  analysis  of 
focal  depths  will  be  performed  for  several  earthquakes  in  the  Kuril/Kamchatka  region. 
Based  on  the  analysis  in  northeast  United  States  and  in  Kuril/Kamchatka  region,  we  can 
begin  to  assess  the  effectiveness  of  these  approaches  for  depth  determination.  We  will 
evaluate  the  recovery  and  identification  of  regional  phases  as  well  as  the  theoretical 
modeling  of  regional  phases. 
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